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ABSTRACT 



In a simplified model we study the hadronization of quark matter in an expanding fireball, in particular 
the approach to a final hadronic composition in equilibrium. Ideal hadron gas equilibrium constrained by 
conservation laws, the fugacity parametrization, as well as linear and non-linear coalescence approaches are 
recognized as different approximations to this in-medium quark fusion process. It is shown that color confinement 
requires a dependence of the hadronization cross section on quark density in the presence of flow (dynamical 
confinement). 



To appear in the Proc. of the Strange Matter 2000 Conference 



CERN-TH/2000-253 
August 2000 



1 Introduction 



Since flavor composition, in particular strangeness enhancement had been suggested as a possible signature 
of quark matter formation in heavy-ion collisions, there is a continuing attention to the question of hadronic 
and quark level equilibrium in these reactions ^|. Recent claims of describing data at AGS as well as at 
CERN SPS by hadronic equilibrium (i.e. using a common volume, temperature and only chemical potentials 
of conserved quantities for the computation of hadron numbers) satisfactorily |j, ||, |j| led to some discrepancy: 
the total reaction time at SPS energy is too short for achieving hadronic equilibrium by hadronic conversion 
processes alone. Experimental data are well described, but poorly understood. After the "go and find it" period 
now we are at the "sit down and think" stage. 

There are also several non-equilibrium approaches to this problem ranging from the use of phenomenological 
parameters on the hadronic level (fugacities) Q through linear and nonlinear (constrained) coalescence mod- 
els Q , and mixed quark- hadron chemistry Q up to very detailed microscopic simulations of the hadronization 
(parton cascades (ToJ, [Tl| VENUS (jj, HIJING jE§, RQMD @ 0, UrQMD @, 0, string model Jl8j, [[§ 



chromodiclectric model [[2l|). The basic physical concepts of these approaches do not include hadronic equilib- 
rium; it might only be reached in course of time. 

The purpose of this talk is to understand the near-equilibrium hadronic result in terms of quark number 
evolution. We consider the simplest possible process on the constituent level: A + B < — ► C. 

In such a process the energy-impulse balance is assumed to be satisfied by the surrounding medium, we 
are interested in the time evolution of the number of final state products Nc(t). The equilibrium number of 
precursors NA,eq and NB,eq - in a quark matter scenario - will be taken as zero. The initial state is assumed 
to be far from this equilibrium, it is oversaturated by A and B with no C particles at all in the beginning. As 
the reaction proceeds there will be C type particles (hadrons) produced. The cross section for such a reaction 
at a first glance can be taken as constant, but temperature and reaction volume must be considered as time 
dependent quantities. We shall point out, however, that the assumption of constant cross sections actually 
contradicts to the confinement principle in quark matter scenarios with transverse flow, so in that case the 
hadronization cross section has to be quark density dependent. This extra dependency, not derivable from 
perturbative or lattice QCD so far, will be referred as "dynamical confinement" in this paper. 



2 Chemistry 

In the chemical approach to hadronization all particle numbers are time dependent. A general reaction can be 
characterized by stochiometric coefficients oq' describing particle number changes in a given reaction (r) in a 
fashion arranged to zero: 

E^ = o, a) 

i 

where Ai symolically denotes one particle of type i. In this arrangement the a^-s can be negative as well as 
positive integers. (For the simple quark fusion process A + B < — ► C one applies ola = &b = 1, olc = — 1.) 

The decisive quantity with respect to approaching chemical equilibrium is the activity of a reaction, 

* w =£^ r) f- (2) 

i 

It describes the (un)balance of gain and loss terms. Rate equations add up contributions from several reaction 
channels, 

# i = £«< r >ijM(l-e* W ). (3) 

(r) 
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Here the factors R^ 1 "' depend on numbers of other particles, cross sections and volume. Due to known (e.g. 
adiabatic) expansion a relation between volume V and temperature T is given. In chemical equilibrium x^ r ' = 
for each reaction, 

E^V* = 0, (4) 

i 

so there are linear constraints for the chemical potentials /li as many as reaction channels. If this is less than the 
number of particle sorts, the equilibrium state is undetermined. If this is greater than the number of particle 
sorts, the reactions cannot be all independent: there are conserved quantities. In this latter case all ^ can be 
expressed by a few numbers corresponding to these conserved quantities. This is the constrained equilibrium. 

Denoting the conserved quantities by Q a and the corresponding charge of the i-th particle sort by qi ia , we 
obtain 

E 9i,aNi = 0. (5) 

i 

Substituting the rate equations (||) we get 

i r 

Since this holds at any instant of an out of equilibrium evolution, the double sum vanishes for arbitrary values 
of reaction activities x^'\ Exchanging the order of summation over particle types i and reaction channels r, and 
taking into account that all factors are positive, we conclude that the considered charges are conserved in 
each reaction: 

£afVa = Qi r) =0. (7) 

i 

In this case we can redefine the chemical potentials as 

M» = E gi.aMa + Mi; (8) 



without changing the activities: 

* (r) = EE^ r) ^+E^ r) §- w 

i a i 

Since the first sum vanishes due to eq. ((?]), chemical equilibrium is equivalent to 

E^ r) f = o. (10) 

i 

For a well-determined or overdetermined system the only equilibrium solution is given by = 0. In this 
(constrained) equilibrium state the original chemical potentials are a linear combination of a few ^ a -s, associated 
to conserved charges, (cf. eq.(^)) 

a 

Whether this chemical equilibrium state can be reached in a finite time depends on the rates and activities 



3 A simple model: A + B < — ► C 



In the case of only one reaction channel, A + B < — ► C, we denote the equilibrium ratio by 

Rcq = 7- • (14) 

iV C,eq 
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Since a a = Q-b = 1 and etc = — 1 for this very reaction, we have Na = Ng = —Nc- From this two conservation 
laws follow, 

N A + N c = N A (0) = N , (13) 

N B + N c = JVb(O) = N . (14) 

For considering an evolution towards equilibrium the specific rate, v = (av) /V, is a useful quantity. It occurs 
in the rate equation 

N c = v (N A N B - R eq N c ) . (15) 
Substituting the conservation laws (|l3|,[l4|) we arrive at 

N c = v ((N - N C )(N - N c ) ~ R eg N c ) . (16) 

This rate is vanishing for the particle numbers Nc = N± with 



N± = R eq /2 + N ± ^R 2 eq /4 + R eq N . (17) 

Here N_ is a stable, N + is an unstable equilibrium point. For an adiabatic expansion R eq = const, for constituent 
hadrons (those with mass rac — tua + Wb). The simple rate equation ( fL6| ) can be solved analytically with the 
result 

: Nc{0) (N c -N + )(Nc-N.)=J t0 Vdt = r - (18) 

Here to is the initial time, when the fireball contains no hadrons, Nc(0) — 0. tb r is the break-up time, 
beyond which no chemical reactions (even no particle collisions) occur. In some cases it is large compared to 
characteristic times and can be approximated by infinity. The solution can be written as 

Nc ^ = n + n - n!-kn. (19) 

with 

K = cxp (- (N+ - 7V_ ) r ) . (20) 

Here N_ is the constrained equilibrium value, achieved for K — 0. There are two alternative forms of this 
relation derived by using 7V + iV_ = N$. The first one is the coalescence form, 

iVc(oo) = K n IZkN- = C ( K ' N oW0)N B (0), (21) 
the second one is the equilibrium form, 

1 - K 
1 - K(N~/N ) 



N c (oo) = 7V_ t ^ /at /at ^ 2 = 7ciVac q . (22) 



In particular Nc (oo) = N- in the K = (infinite rate) case. On the other hand in the case of confinement 
R eq = 0, from which N + = 7V_ = A^ and K = 1 follow. 

For (N + — N-) small (0(R eq )), but non-zero the precursors are already hadronic, but suppressed in the 
equilibrium state. Then we expect, that in the R eq — ► limit the precursors A and B are no more present 
in the final mixture. This is, however, insufficient for the quark matter hadronization scenario, because the 
above limit does not automatically lead to an absence of A and B particles in the final state, if the expansion 
is dramatic enough. We describe this phenomenon in the followings. 

Let us expand K in terms of (N + — iV_) which is in the same order of magnitude as R cq — > 0. We get 

K = exp(-(N + - N-)r) « 1 - r(N+ - N-), (23) 

and the final hadron number becomes 

N C (^) = N -^-, (24) 
1 + rNn 
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leaving us with a non- vanishing precursor number 

N A (oo) =N B (oo) = N 1 . (25) 
1 + rN 

Eq. (|i|) has different readings. For small r it is close to the result of linear coalescence models 

N c (oo) = rNl (26) 

(Note that N A (0) = N B (0) = No.) ALCOR uses extra modification factors b A and 6b in the coalescence rate 

N c {oo) =rb A b B Nl (27) 

To start with it seems to be rather similar to the coalescence picture. Due to conservation of quark number 
during an - supposedly rapid - hadronization, however, ALCOR obtains 

N c (oo) = rb A b B N^ = —^Nt (28) 
1 + rN 



and determines the factors bi accordingly 



Vl + rN 



It leads exactly to the constrained equilibrium result in the infinite integrated specific rate limit r Nq — ► oo (then 
b A = b B w 1/Vr^o): 

N c {oo) = N , N A (oo) = N B (oo) = 0. (30) 

In this limit the precursor numbers approach zero, and the above discussion can be applied for quark matter 
hadronization. Or, reversing the argument, we conclude that for quark matter hadronization r has to diverge 
in order to agree with confinement in the end stage. (This fact is nonetheless hidden in the thermodynamical 
limit with infinite quark number No — * oo.) 

Finally the chemistry approach gives the interpretation of the fugacity due to 

N c (oc) = iV 2 y^j^ = ^cNo- (31) 

In fact the relationship between these approaches is more pronounced by inspecting chemical potentials instead 
of final numbers. In the discussed simple case the chemical potential of the (final state) hadron is given by 

f=log^. (32) 

J Nc,eq 

Substituting our result for Nq (oo) it splits to contributions with different physical meanings: 

— = log h log = — ! ■ h log 7c- (33) 

T N C: eq N- T 8 '° V ' 

The conserved quantities (charges) carried by the hadron C are those involved in the identities respected by 
the rate equations, N A (t) + N c (t) = N A (0) = N and N B (t) + N c {t) = N B (0) = N . In this case q c ,i = 1 and 
qc.2 = lj respectively. The fugacity parameter (as used by Rafelski et al.) is the remaining part of the formula, 
in equilibrium it approaches one: 7c = 1. 

We note here that the ALCOR approach (constrained coalescence) leads to the equilibrium result only if 
a single reaction channel is considered, like in the case discussed in the present paper. For several, competing 
channels the relative hadron ratios are different for ALCOR and constrained equilibrium. This best can be seen 
by inspecting the corresponding chemical potentials to the pre resonance decay hadron numbers as a function 
of strangeness in the ALCOR simulation. Slight deviations from the perfect straight lines, which would mark 
the constrained equilibrium in the ideal case, can be observed in particular for multiply strange anti-baryons 
(cf. Fig.l). 



4 




Figure 1 : Equivalent chemical potentials of composite hadrons before resonance decay obtained in ALCOR as 
a function of the strangeness. Linearity corresponds to hadronic equilibrium. 



4 Expansion effects 



We continue by obtaining the integrated specific rate, r, for different expansion scenarios. This quantity depends 
not only on cross section(s) - first assumed to be constant and later density dependent - but also on the expansion 
and cooling of the reaction zone. We consider here two cases: one dimensional expansion with a rather relativistic 
equation of state and three dimensional expansion of a massive, non-relativistic ideal Boltzmann gas. The first 
case is intended to be characteristic for RHIC and LHC, the latter rather for SPS and AGS energies. 

During a one-dimensional expansion the volume grows linearly in time, V = irR^t, and the (thermal) average 
of relative velocities is near to the light speed (v) = 1. In this case we obtain 

a dt (avo) t br 

^7 = ^ olog V (34) 

The time-integrated specific rate grows with the break-up time t br monotonically, it can eventually be infinite, 
if there is no thermal break-up of the expanding system. In that case 7c = 1, equilibrium is achieved. For a 
three-dimensional expansion one gets V = Vo(t/to) 3 . Furthermore a massive, non-relativistic ideal gas satisfies 
VT 3 / 2 = V Tq /2 during an adiabatic expansion. The average relative velocity is temperature dependent like 
(v) ~ y/T/m ~ 1/t, with some reduced mass of the fusing pair m. This together leads to the estimate 

Evaluating the integral we arrive at 

1 (av ) 



3 V 




to i- r • ( 36 ) 



It is important to note that r is finite even without break-up, in the t^ — > oo limit. This means that the 
non-relativistic, 3-dimensional spherical expansion scenario with constant hadronization cross section cannot 
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comply with quark confinement. This result is valid in the presence of any, however mild, transverse flow, since 
the one-dimensional expansion scenario showed a logarithmic dependence on the break-up time. 



Let us now consider a quark density dependent hadronization cross section which is at least inversely pro- 
portional to the density of quarks (n ~ \jy ^ i -3 ). Exactly this assumption was made in transchemistry B. 



Introducing the initial collison time, 
with density po = 2Nq/Vq we arrive at 



(<™o)Po 



aVo) ' (38) 



V 2N t c 

The quantity rN describing the deviation from equilibrium for hadrons then simplifies to 

WV = ^log^ (39) 



2t c to 



3 X 



rJVo = #- l-(^) I (40) 



for one-dimensional relativistic flow and 



6i c V \tb, 

for three-dimesional non- relativistic flow of massive constituents. From these formulae one realizes that for 
ttr = to, i-c in case of immediate break-up the chemical equilibrium cannot be approached, no hadron C 
will be formed and 7c remains zero. The opposite limit, t\> r — > 00 is less uniform: while in the relativistic 
one-dimensional scenario equilibrium is always reached if enough time is given (7(7(00) = 1), this is not so for a 
stronger, more spherical expansion or in the presence of a transverse expansion. In the case of three-dimensional 
spherical expansion 7c tends to a finite value, 

7c(oo) = 1 , (41) 
1 + 6t c /t 

which is in general less than one. It means that chemical equilibrium cannot be really reached even in infinite 
time. (An analysis of proton - deuteron mixture with a similar conclusion was given in Ref.[p2|.) For early 
break-up, ibr — to <C to, both scenarios lead to 

7c « « 1. (42) 



In quark matter hadronization scenarios eq.(|4l|) contradicts to the quark confinement principle. In order to 
reach 7c (00) = 1 a dynamical confinement mechanism, causing the quark fusion cross section in medium to 
be (at least) inversely proportional to the quark density, has to be taken into account besides the equilibrium 
suppression of quark number (referred in this paper as "static confinement"). 

Let us assume that A and B are colored objects, and the in-medium quark fusion cross section scales as 

V 2N x 1+e 



'-"{knits:) (43) 

with some positive e. The solution of the rate equation in the R eq = case (both static and dynamic confine- 
ment) becomes 

N c = N (l - (1 - S) 1 /^ , (44) 

with 

1 (45) 




G 



Equilibration 



RHIC dynamic 



SPS dynamic 



0.8 



O 0.6 




1 1 1 1 1 1 1 1 1 1 

123456789 10 

>br 'to 

Figure 2: Degree of the chemical equilibrium reached for different break-up times ibrAo at initial inter-collison 
times t c /to = 1/3 for SPS (3-dimensional) and RHIC (1-dimensional) expansion scenarios with static and 
dynamic confinement (hadronization cross section is inversely proportional to the color density). 



in the 3-dimcnsional flow (SPS) scenario. For e = — 1 (without dynamical confinement) one gets back the result 
discussed so far. In the limiting case e = (minimal dynamical confinement) one obtains 

Nc = N (l- (7)^) ■ ( 46 ) 

For a 1-dimensional flow (RHIC) scenario we arrive at 

N c =N (l- e-t*-^)/ 24 ^ (47) 

in the e = case. 

Fig. 2 plots the dependence of the non-equilibrium ratio Nc/Nq on the scaled break-up time t& r /to f° r these 
scenarii. 



5 Conclusion 



In conclusion 1) the results of the statistical model can be interpreted not only by assuming hadronic equilib- 
rium, but as well by quark - hadron mixture chemistry with quarks eventually suppressed due to confinement. 
This interpretation avoids the collision time problem: the hadrons need not transmutate into each other directly. 
2) The ALCOR and constrained hadronic equilibrium leads to the same result in the simple scenario presented 
here due to the constraints on conserved quantities. Quark coalescence, which takes conservation into account, 
closely mimics hadronic statistics in chemical equilibrium. This feature does not generalize for several com- 
peting reaction channels, but the overall picture is still similar. 3) 7c ~ 0.3 is typical for spherical expansion 
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(CERN SPS) without considering dynamical confinement, even for an arbitrary long hadronization process. 
The experimental evidence, reaching at least 70-80% of equilibrium ratios for almost all hadrons in a short 
time in the order of several fm/c, excludes this possibility. It also underlines the quark matter hadronization 
scenario with dynamical confinement effect, as transchemistry utilizes it. 4) A purely longitudinal expansion 
scenario, perhaps typical for RHIC and LHC energies, in principle allows for reaching the equilibrium even 
with static confinement only, but the process is rather slow and an early break-up may prevent some quarks 
from hadronizing. Dynamical confinement has to be in work for near one-dimensional expansion scenarios with 
a minor transverse flow, too. 5) We pointed out that total elimination of free quarks in an expanding quark 
matter scenario in the presence of any transverse flow requires a quark density dependence of the hadronization 
cross section, which is at least 1/n. We call this effect "dynamical" confinement, it evidently cannot be of 
perturbative origin. 

One is tempted to speculate about the dynamical confinement, looking for a cause of the density dependence 
of the hadronization cross section. Since this dependence makes the cross section stronger and even diverging 
with diminishing quark density, it cannot have perturbative origin. In our view it has to do with the formation of 
strings or color ropes, which become longer and hence more efficiently bind quarks into hadronic (color neutral) 
clusters as the color charge density drops. 
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